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ABSTRACT 

We present cross sections and reaction rates for merging to occur during encounters 
of equal-mass spherical galaxies. As an application, we determine the rate of galaxy 
merging in clusters of galaxies. We present results for two types of Plummer models (a 
full and a truncated one), two King models and the Hernquist model. Cross sections 
are determined on the basis of a large number (~ 500) of simulations of galaxy encoun- 
ters, using the 10-Gigaflops GRAPE 3A special-purpose computer. We characterize the 
overall merger rate of galaxies in a galaxy cluster by a single number, derived from our 
cross sections by an integration over galaxy encounter velocities in the limit of a constant 
density in velocity space. For small clusters, where the cluster velocity dispersion may 
not significantly exceed the internal velocity dispersion of the individual galaxies, this 
constant-density approximation may not be valid. For those cases, we present separate 
results, based on integrations of our cross sections over Maxwellian velocity distributions. 
Finally, tidal effects from the cluster potential as well as from neighboring galaxies may 
prevent a barely bound galaxy pair from spiraling in after their first encounter. We give 
a quantitative estimate of the resulting reduction in the actual merger rate, due to these 
tidal interactions. 
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1. Introduction 

Detailed simulations of galaxy encounters have be- 
come increasingly sophisticated during the last twenty- 
five years, due to significant improvement both in 
computer hardware and in algorithms used (for a re- 



cent review, see Barnes and Hernquist 1992). Most of 
these simulations attempt to explain features of spe- 
cific encounters, often in an attempt to reproduce par- 
ticular observations of interacting galaxies or merger 
remnants. In contrast, encounters between simpler 
galaxy models in order to obtain more general statis- 
tical results, have received less attention. 

The present paper attempts to address this lat- 
ter problem, by limiting ourselves to simulations of 
equal-mass galaxy models of various types, in order 
to determine merger cross sections and reaction rates. 
Our results can be readily applied to study the evo- 
lution of galaxy clusters, if we use our models to rep- 
resent dark matter halos around galaxies in the limit 
that these halos can be considered to be spherical. 
In practice, moderate deviations from halo sphericity 
will not greatly afi^ect our results. 

The merger probability in an encounter of two 
galaxies is enhanced significantly when the encounter 
takes place at relatively low speed. For a given clus- 
ter of galaxies, with a given population of specified 
galaxy halos, the net merger rate can be determined 
by an integration of the usual nav factor (density x 
cross section x velocity) over a Maxwellian or low- 
ered Maxwellian velocity distribution. This proce- 
dure is sketched qualitatively in Fig. lb, and can be 
compared to that used in stellar evolution calcula- 
tions (Fig. la) where nuclear reaction rates are de- 
termined from the encounters of individual nuclei in 
the cores of stars. In the latter case, cross sections 
rise quickly with increasing encounter speed, leaving 
a small window for effective encounters well into the 
high-energy tail of the Maxwellian distribution. Our 
case of galaxy mergers, in contrast, is dependent on 
low-speed encounters, and determined by a relatively 



small window at the low-energy end of the Maxwellian 
distribution. 

Fig. lb suggests that we can make a good start 
by representing the velocity distribution of the galax- 
ies in the cluster as having a constant density in ve- 
locity space, i. e. f{v) cx in three dimensions. 
This constant density can be obtained, for any given 
Maxwellian distribution, as the density of galaxies 
around the peak of the folded merger occurrence fav 
in Fig. lb. In §2 we describe some technical details of 
our calculations. In §3 we present the cross sections, 
from which we determine reaction rates in §4. In §5, 
we determine the corrections to the reaction rates, 
due to tidal effects. As an application, we present the 
rate of galaxy merging in clusters of galaxies in §6, 
and compare that with results from previous papers. 
§7 sums up. 

2. Simulations 

2.1. Units 

We choose our physical units of mass, length, and 
time by requiring that AI = = G = I, where M is 
the mass on a single galaxy, is the virial radius of 



a galaxy and G is the gravitational constant (Heggie 



and Mathieu 1986). The virial radius is a measure of 
the size of a galaxy, defined so as to be able to express 
the potential energy of the galaxy to be equal to 

GA'P 



Epot — 



2r„ 
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In practice, the virial radius is not very different from 
the half-mass radius r^. Typically, 

Th =i 0.8r„, (2) 

with a coefficient generally in the range 0.76 - 0.98 



(Spitzcr 1987, §1.2a). In our standard case of a Plum- 
mer model we have 
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and for a constant-density sphere 
3 X 22/3 



rh = 



-Ty = 0.95ri,. 



(5) 



The main reason to take the virial radius, rather 
than the more usual half-mass radius, as our choice 
of units is a physical one: when we vary the internal 
structure of our galaxies, it is more meaningful to 
keep the total energy £^ of a galaxy constant than its 
half-mass radius. In our units. 



E = Epot/^ = -1/4. 



(6) 



In terms of the particle positions, the virial radius 
is the harmonic mean particle separation 
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averaged over all particle pairs i ^ j. 

2.2. Galaxy models and initial conditions 

For the initial spherical galaxy model, we adopted 
five different choices: two Plummer models with dif- 
ferent cutoff radii Tc (22.8 and 4), two King mod- 
els with different values for the central potential Wc 



(1 and 7), and a Hernquist model ( Hcrnquist 199C ) 
with Tc — 20. We will indicate the Plummer model 
with = 4 by 'truncated Plummer' and the Plmnmcr 
model with = 22.8 simply by 'Plummer'. 

These models provide a wide range of central 
condensation and halo distribution. They all have 
an isotropic velocity distribution since, according to 
Aguilar and White (1985), the difference in the ve- 
locity distribution is less important in determining 
merger cross sections than the difference in the den- 
sity distribution. For all calculations, initial galaxies 
are modeled by 2048 equal-mass particles (except for 
several test calculations, see §3.4). Thus the total 
number of particles per run is 4096. For Plummer 
models and King models, we constructed random re- 
alizations using an algorithm similar to the one de- 
scribed in Aarseth, Henon and Wielen (1974). For 



Hernquist models, we adopted a kind of "quiet start", 
in which particle i is initially placed in an arbitrary 
position on a sphere with a radius corresponding to 
the Lagrangian mass (i — 0.5)/7V. The velocity of each 
particle was then chosen randomly from the velocity 
distribution function at that radius. 

Since the initial galaxy models are spherical, the 
relative orbit of the two given galaxies is determined 
by only two parameters specifying the incoming branch 
of the hyperbolic Kepler orbit, the impact parameter 
p and the relative velocity at infinity v. In this study, 
we found it more convenient to specify the initial orbit 
by V and the pericenter distance of the extrapolated 
unperturbed hyperbolic orbit, in the limit that two 
galaxies would have been point masses. The reason 
to use Tp rather than p is related to the measurements 
of energy change during relatively wide encounters: 
keeping Tp constant allows us to let v smoothly go to 
zero in the parabolic limit, where p — s- cx). 

Fig. 2 shows the initial conditions used to deter- 
mine the merging criterion for each models. The num- 
ber of runs for one model is 60-140. The following val- 
ues for Tp were used: 0, 0.0125, 0.05, 0.2, 0.4, 0.6, 0.8, 

1.2, 1.4, 1.6, 1.8, 2.0, 2.4, 2.8, 3.2, 3.6, 4.0, 4.8, 5.6, 
6.4, 8.0. For the initial separation between the two 
approaching galaxies, we choose a value of roughly 
10 max(rp, 2r?i), where is the half mass radius of 
each of the two identical galaxies. 

2.3. Hardware and Numerical Integration Method 

We have used a GRAPE-3A ( |Okumura 1993| ), 



a special-purpose computer for coUisionless A^-body 
simulation, for all simulations. The GRAPE-3A per- 
forms the calculation of the gravitational interaction 
between particles. It is used with a host computer 
which is a UNIX-based workstation. All calculations 
except for the force calculation are performed on the 
host computer. One GRAPE-3A board hosts eight 
custom LSI chips, each of which calculates the grav- 
itational interaction between particles in the speed 
of about 0.75 Gflops. Thus the peak speed of one 
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board is 6 Gflops. For most runs, we have used two 
GRAPE-3A boards in parallel. 

For the force calculation, we have used a simple 
direct-summation algorithm. Although it is perfectly 
possible to use the O(A^logA^) tree algorithm on the 



GRAPE systems (Makino 1991), it is not practical to 



do so for small particle numbers. The minimum num- 
ber of particles for which the tree code is faster than 
the direct summation algorithm is 1 ^ 4 x lO'^, de- 
pending on the speed of the host computer. Since we 
typically used only 4096 particles, direct summation 
was much faster than the tree code. 

We used the time-centered leap-frog integrator. 
The time step was 1/128 for all runs. The size of the 
softening parameter was e = 1/32. These parameters 
gave sufficient accuracy, resulting in a relative energy 
error well below 0.1 % at the end of each calculation). 

2.4. Simulation Termination Procedure and 
Data Reduction 

For each choice of galaxy model, and for each 
choice of pericenter distance r^, we have carried out 
a series of galaxy encounter simulations at different 
values for the encounter velocity v (measured at in- 
finity). For each of these simulations, we determine 
the energy of the relative orbit of the two galaxies, 
well after the first encounter, using the simple energy 
criterion: 



E 



GM1M2 



1 



-flV 
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where Mi and M2 are the masses of two galaxies, r is 
the distance between two galaxies, G is the gravita- 
tional constant, ji = M1M2/ (Mi -I-M2) is the reduced 
mass, V is the relative velocity. 

In Fig. 3, we plot the energy E{t) as a function of 
time t, for the case of Hernquist model initial condi- 
tions, with Tp — 0.0125 and v — 1.2. It is clear that 
the asymptotic value of the energy is reached soon 
after the encounter. However, to be on the safe side, 
we have extended all runs to reach a final separation 
comparable to the initial separation, where possible 
(for a bound final state, we have taken the minimum 



of this separation and the apocenter of the elliptic 
orbit). In this particular example, it takes far more 
time to reach the initial separation again, because the 
outgoing orbit has far less positive energy and there- 
fore the galaxies recede more slowly than they came 
in. 

Having determined the asymptotic energy values 
E{rp,v), we select the pair of v values between which 
the energy changes sign. We then use linear interpo- 
lation to determine the critical point rp{v) at which 
the outgoing orbit would have been just parabolic, 
separating the regions of merging and escape. Fig. 2 
indicates this process, for the case of Plummer mod- 
els. Each circle represents a galaxy encounter run, 
and the full line connects the critical rp{v) points. 

Fig. 4 shows the same results, but translated from 
Tp values to p values. The impact parameter p is re- 
lated to the pericenter distance Tp through the gravi- 
tational focusing relation (see Eq. |l^ below) . 

As a technical note, here is a brief description of 
the actual implementation of the energy determina- 
tion, Eq. (^, in which appear the mass, position and 
velocity of the galaxies after the encounter. They are 
calculated by the following procedure: 

1) Make a crude guess for the center of mass of each 

galaxies. Here, we use the center of mass of 
particles that are initially in one galaxy as the 
guess for the center of mass of the galaxy at that 
time. 

2) For each particle, determine which galaxy it be- 

longs to. We calculate estimated binding energy 
given by 



GMj 1 2 



(9) 



where rjj is the distance between particle i and 
the estimated center of mass of the galaxy j, Vij 
is the relative velocity between them, and Rq is 
the length scale parameter that represents the 
depth of the potential well of galaxies. We use 
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Rq = 0.6, which gives a fairly accurate estimate 
for the potential for Plummer models. We then 
make our initial guess as to which galaxy parti- 
cle i mostly likely might be bound, by assigning 
it to galaxy 1 if En < Ei2, or to galaxy 2 if 
En > Ei2. 

3) For each galaxy, determine which of the particles 

that are labeled as belonging to it are actually 
bound to the galaxy. If the binding energy of 
a particle relative to the galaxy is positive, we 
regard it as an unbounded escaper. 

4) Repeat step (3) until membership converges. 

Most of the data reduction was performed using 
the NEMO software package, in the version provided 
by Peter Teuben (1995). 

3. Cross Sections for Galaxy Merging 
3.1. Results and Scaling Relations 

Fig. 5 shows the critical velocity for merging Vcrit 
as a function of the impact parameter p, for all galaxy 
models used. The merging criterion determined ex- 
perimentally here shows only a weak dependence on 
galaxy model: for the whole range of impact parame- 
ter studied, the difference in Vcrit among the different 
models is less than 20%. For the Plummer model, we 
find a good fit with Vcrit oc p^^-''^ for large p, as plot- 
ted for comparison in Fig. 5. For the King model with 
Wc = 1 the index is slightly larger and for Wc = 7 the 
index is slightly smaller. The slope for the Hernquist 
model is somewhat more small. 

At first sight, these results seem to be in good 
agreement with the distant tidal impulsive approxi- 
mation, which indeed predicts a relationship Vcrit oc 



P 



-0.75 



for large p. The tidal impulsive approximation 
has been used by many researchers to obtain analytic 
approximations for changes in mass and energy of 



galaxies involved in relatively weak encounters ( Rich- 



tone 19751 , [White 1979| , [Dekel, et al. 1980| , [Richstonc 



of energy of the relative orbital motion of the two 
galaxies scales as 



\rivp 



(10) 



where rp and Vp is the relative distance and relative 
velocity at closest approach. For the critical velocity, 
the energy loss is equal to the kinetic energy of the 
orbital motion at the infinity: 

1 



(11) 



where p — M/2 is the reduced mass. If we apply the 
point mass approximation for the relative orbit, the 
relations between the orbital elements are given by 
conservation of angular momentum and energy as 



pv = rpVp, 



-pv 



1 2 GM^ 
r2 , 



(12) 



Here, we assume that pv'^ << GK'P /vp. Under this 
assumption, the last equality is reduced to 



/4GM 



(13) 



The resulting pericenter distance is then given as 



? 2 

p^v^ 



(14) 



^ 4GM' 

and the incoming velocity can now be expressed in 
terms of the impact parameter as 



(15) 



md Malumuth 1983). In this approximation, the loss 



seemingly in good agreement with the data shown in 
Figure 6. 

3.2. The Masking Tendency of Gravitational 
Focusing 

This good agreement is, however, fortuitous and 
does not reflect the physical reality. The discrepancies 
can be seen clearly in Fig. 6, which shows the critical 
velocity Vm as a function of pericenter distance Vp. 
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If the tidal approximation were really valid, the re- 
sults for all models would show the same asymptotic 
behavior oi v oc Vp^^'^ . The results for the Plum- 
mer models and the concentrated King model is close 
to the theoretical line. However, the results for the 
Hernquist model and the shallow King model deviate 
significantly. 

Indeed, there was no reason to expect an impulsive 
tidal approximation to give us guidance in estimating 
the critical velocity for mergers to take place. As is 
clear from Fig. 6, mergers at large Tp take place for 
nearly parabolic orbits, where the duration of rela- 
tively close encounters around pericenter passage is 
drawn out to a total time exceeding that of the half- 
mass crossing time. There is no reason to expect, un- 
der these circumstances, that any type of impulsive 
approximation would be valid. 

Why, then, did Fig. 5 show such a good correspon- 
dence with the line Vcrit oc The reason can be 
found in a conspiracy of gravitational focusing effects 
together with the real merger criterion. This can be 
seen as follows. 

Suppose that we start with a relation for the energy 
loss in the form of a power law 

A£;cxrp". (16) 

For the critical velocity at infinity, we then have 

i)cxr;"/^ (17) 
From Eq. (|l^, we obtain 

wcxp^. (18) 

The slope of this relation between p and v lies in the 
range {—1,0}, for any positive value of a. This im- 
plies that even large differences in the relations be- 
tween nearest approach and energy loss are not very 
well visible in Fig. 5. 

From a practical point of view, this reasoning has 
an interesting consequence: Eq. ( p^ ) implies that 
the merging cross section is relatively insensitive to 



the fact that the amount of energy loss varies widely 
among different models, when measured at the same 
pericenter distance. 

From a theoretical point of view, we already saw 
that the impulsive tidal approximation is not reason- 
able. More specifically, the time scale of the encounter 
is proportional to r^'^. If this encounter time scale sig- 
nificantly exceeds the orbital time scale of a particle 
in the galaxy, the binding energy of that particle be- 
comes an adiabatic invariant. As a result, the energy 
change decays exponentially as a function of increas- 
ing Tp. If a galaxy has a finite radius r, encounters 
with Tp 3> r will result in only a very small energy 
change, and v{rp) would thus decrease exponentially. 

3.3. Theoretical Scaling Arguments 

A simple alternative to the impulsive approach 
could start with a determination of the amount of 
significant particle overlap during the encounter. A 
quick inspection of Figs. 5 and 6 already shows that 
this may be a reasonable approach, when we real- 
ize the exceptional nature of the Hernquist and King 
(Wc = 1) models. These are the two models that de- 
viate most in both figures, and indeed, these are the 
two models that deviate most in the position of the 
radius that encloses 95% of the mass in each model, 
with the King {Wc — 1) model having an unusually 
sharp cut-off and the the Hernquist model having an 
unusually large amount of matter far from the core 
region (see also Fig. 12). 

If a galaxy has an extended halo with a power-low 
density profile, the particles in the region outside the 
radius rp have orbital time scales comparable to, or 
longer than, the time scale of the encounter. As a 
result, those particles are strongly perturbed, and we 
can still expect the impulsive approximation to give 
a reasonable estimate, at least in order of magnitude. 
Thus the energy gain of particles outside the radius 
Tp can be estimated from the typical velocity kick 
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they receive as 

f 2GM \ ^ 

AEout ~ [M-A/(rp)] {^vf = [M~M{rp)] 

V TpVp J 

(19) 

where M{r) is the mass of a single galaxy enclosed 
inside the radius r (the total mass A/ = 1 initially, 
in our units). The velocity Vp at pericenter passage is 
related to Vp in the parabolic approximation by « 
GM/rp. This gives: 



G[M - M{rp)]M 



(20) 



Note that this energy change is equal to the binding 
energy of the mass outside the radius Vp. In other 
words, in our approximation all particles are divided 
into two groups, an inner and an outer group. The 
inner particles are considered to be left undisturbed, 
while the outer ones are significantly perturbed dur- 
ing the encounter. 

For galaxy models with a halo described by a 
power-low density pd — , we have [M ~ Af (r)] cx 



AEout oc r: 



V (X r. 



2-13 
P ' 
1-/3/2 



(21) 



For the Hernquist model, /3 — 4. Therefore the power 
index in the latter relation is —1. For the Plummer 
model, (3 — 5 and the index becomes —3/2. 

Fig. 7 shows the critical velocity obtained by A''- 
body simulations and that estimated using Eq. (|20|). 
For both Plummer and Hernquist models, the theo- 
retical estimate and the numerical result show a rea- 
sonable agreement, well within a factor two for Vp val- 
ues in the range {1,20}. This is a very satisfactory 
result, given the simple nature of our approximations. 

3.4. Cross Sections 

From the Pm{v) results, plotted in Figs. 4 and 5, it 
is straightforward to determine the cross section a{v) 
for mergers to occur: 



The results are plotted in Fig. 8. 

A more useful way to display these cross section 
is by multiplying them with a factor v'^. The first 
two factors of v reflect the three-dimensional nature 
of a Maxwellian velocity distribution, with a velocity 
space volume factor of v^dv (see Eq. |2^), while the 
third factor indicates the fact that, for a fixed target 
size, the rate of merging encounters is proportional to 
the relative velocity of the galaxies involved. 

The results, in the form of v^a{v), are plotted in 
Fig. 9. There is a slight gap at the left-hand side 
of the curves, where it would have been too time 
consuming to determine the critical merger velocities 
through numerical simulations. Fortunately, we can 
use the arguments developed in §3.3 to extrapolate 
our numerical results leftwards of where the curves 
end. For galaxy models with a halo described by a 
power-low density pd = r^^, we find from Eqs. (14) 
and (|2ll) that 



p (X V 



This leads directly to 



(,3-l)/(2-/3) 



(T«3cxpV(X«(^-4)/(/'-2). 



(23) 



(24) 



For a Hernquist model, with /? = 4, we see that this 
quantity tends to a constant value for vanishing v, 
while for a Plummer model this quantity tends to 
scale as v^^^ for low v values. In the rate determina- 
tions of §4, we will use these analytic extrapolations 
to augment the numerical data. 



3.5. 



Error Discussion 



In our simulations, the total number of particles 
per galaxy has been typically N = 2 x lO'^. Depending 
on the exact definition of half-mass relaxation time 
thr and half-mass crossing time t^c, we can get some- 
what different relationships between these two quan- 
tities. Let us take as a reasonably accurate choice 
the expression thr/thc — O.liV/ ln(0.4A^) (see 3pitzei 



(t{v) — ■np^ 



(22) 



1987 for the factor 0.4). For iV = 2 x 10'\ this gives 
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us thr/thc — 30. However, we have used a soften- 
ing length of 1/32. With a strong-deflection distance 
of order the added softening increases the re- 

laxation time by a factor log TV/ log 32 ~ 2. Our 
galaxies can thus be expected to show significant re- 
laxation after 60 half-mass crossing times. In our 
units, the internal three-dimensional velocity disper- 
sion of a galaxy is l/\/2, and the time to cross the 
system, starting at the half-mass radius, is therefore 
~ 2/^/2 ~ 3. Thus, we can expect significant relax- 
ation to occur after 200 time units. 

For the simulation in fig. 3, relaxation efl[ects are 
not expected to be very important, given a total du- 
ration of the simulation of 100 time units. However, 
for larger impact parameters, the duration of a typical 
encounter simulation can be significantly longer, and 
relaxation effects may begin to influence the measure- 
ments of energy dissipation. This in turn will effect 
the determination of the border line for mergers to 
take place at high values for the impact parameter. 

We have run various tests in order to determine the 
systematic errors that can occur for large Vp values. 
For example, for the King models with Wc = 1, we 
have rerun a number of experiments, using a total 
particle number per galaxy in the range 512 < N < 
16384, as compared to our standard value N = 2048. 
For relatively small pericenter values, such as Vp = 0.1 
and Tp = 3.2, we did not find any noticeable change 
for the final merger boundary value Vm{rp)- In fact, 
even downgrading our runs by a factor four in total 
particle number did not have a significant effect on 
the outcome. 

However, for Vp = 6.4, we found that the mea- 
sured Vmifp) steadily decreased for increasing parti- 
cle numbers, until leveling off for a value of N four 
times larger than originally used. As a consequence, 
the reduction in slope in the lower right corner in the 
line labeled King (Wc = I) in Fig. 6 is not correct. 
In the limit of ^ oo, the line will continue to de- 
scend with a near-constant slope, leading to Vm values 



around rp = 8 that are roughly half of what we have 

obtained here for our standard value of N = 2048 
particles per galaxy. However, these systematic errors 
do not contribute significantly to the overall reaction 
rates, as we can see in Fig. 9. 

4. Merger Rates 

Given the cross section determined in the previ- 
ous section, we are now in a position to determine 
the rate at which mergers occur, by averaging over a 
distribution of velocities. 

The distribution function for the velocity of the 
galaxies inside a cluster is given in a Maxwellian ap- 
proximation by 

/iW=(^)'^^e-Ve-V2-^, (25) 

where (Je is the one-dimensional velocity dispersion 
for the cluster. 

The distribution of encounter velocities is then 
given by 

Mv) = 2-'n-'/'a-'v^e--'/^-", (26) 

where in both cases we use a normalization condition 
of unity after integration over velocities: 

/•OO POO 

/ fi{v)dv= / f2{v)dv = l. (27) 
Jo Jo 

Let us introduce the ratio x of the one-dimensional 

velocity dispersion of the cluster and the one-dimensional 

internal velocity distribution of the stars in each 

galaxy: 

X= — = \/6£7e, (28) 

since crj = l/\/6 in our units in which the specific 
kinetic energy ^af = 1/4. We thus have 

h{v) = 2-\-^/^a-\^e-i^"/^" , (29) 

The number of merging events per unit time per 
unit volume is given by 

/•OO 

R{n,x) = n^ vf2{v)a{v)dv, (30) 

Jv=0 
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where n is the number density of the galaxies and 

a{v) is the cross section of merging. 
We can also write this result as 

R = n'^2-\-^/^a-^R{x), (31) 

where R{x) is the non-dimensional merging rate de- 
fined as 

R{x)= r"\^a{v)e-i'''/'''dv, (32) 
Jo 

where the upper limit of integration is reached at Vcrit, 
since the cross section a{v) = for v > Vcrit- 

Note that in our units it is not immediately clear 
that this last expression is dimcnsionlcss. If we re- 
mind ourselves of the fact that our internal one- 
dimensional velocity dispersion ai = l/\/6 and the 
virial radius of an individual galaxy = 1, we can 
rewrite the expression for R{n, x) as 

R = n'2-'7T-'/^rl^^^R{x) = ^l^n^rlaiR{x), 

<J^ -^TT X^ 

(33) 

from which it is clear that R{x) is dimensionless, since 
all physical factors (density squared x cross section 
X velocity) have been displayed here explicitly. 

In rich clusters, the velocity dispersion ae in the 
cluster is several times larger than the internal veloc- 
ity dispersion cr^ in each of the galaxies. Even for less 
rich clusters, ae > <7i, typically. 

Thus, for most rich clusters, the merging rate cal- 
culated with constant f{v) gives an error less than 
a few percent. For that regime, the dimensionless 
merging rate is given by 

Roo = lim R{x) = / v^a{v)dv. (34) 

Note that _Roo docs not depend on the velocity 
distribution of galaxies within the cluster. The only 
dependency on the external velocity dispersion cTj is 
through the factor in the expression for R itself. 
This is simply the dilution factor in velocity space: 
increasing the external velocity dispersion decreases 



the density of galaxies within the region in which en- 
counters can lead to merging. 

Table 1 shows the nondimensional merging rate 
i?oo for all models. For non-zero values of x, we can 
approximate R{x) for most models by the following 
fitting formula 

while for the Hernquist model a better approximation 
is given by the fitting formula 

The physical reason behind the fact that the Hern- 
quist model shows a merging rate that is 15% larger 
than the Plummer model lies in the more extended 
mass distribution of the former, as is illustrated in 
Fig. 12. 

5. Tidal effects 

In our discussion of the merger cross section and 
the merger rate, we have so far regarded two galax- 
ies as actually merging if the orbital binding energy 
of the two galaxies after the first encounter is nega- 
tive. Simple and straightforward as this assumption 
is, it may not be realistic in some extreme cases, as 
for example in rich clusters. In such cases, even if 
two galaxies have become formally bound after their 
first encounter, they may still be disrupted under the 
influence of either the tidal forces from the cluster as 
a whole, or the tidal forces from individual nearby 
galaxies. In the following, we evaluate the magnitude 
and relevance of both of these effects. 

5.1. Tidal effects from the cluster as a v^rhole 

Consider a cluster consisting of N galaxies. The 
tidal force from the cluster would dissolve a pair of 
galaxy if the separation of the pair is larger than 
fR[ND)~^/^, where R is the radius of the cluster 
and D is the ratio of the total mass of the cluster 
to the mass associated with the individual galaxies. 
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The coefficient in tliis expression would be / = 1/3 
in the Umit in which the cluster mass is all concen- 
trated in the center (Spitzer, 1987, eq. 5-5). For more 
realistic mass distributions, the tidal force on a typi- 
cal merging pair is considerably less, and therefore / 
is significantly larger. A reasonable approximation is 
therefore to simply take / = 1. 

Thus, we discard a candidate merger pair if their 
first apocenter distance exceeds 



5.2. 



rt..giobai = R{ND)-^'\ (37) 
Tidal effects from individual nearby galeLX- 



To determine the tidal effects of neighboring galax- 
ies, acting as perturbers, requires a somewhat more 
complicated calculation. For a pair of galaxies in the 
process of merging to become unbound by the en- 
counter with a third galaxy, this latter galaxy should 
pass within a distance that is sufficiently close, during 
a time that is sufficiently quick, i.e. before the first 
fall-back. Since most orbits after the first fly-by are 
highly eccentric, the fall-back time can be approxi- 
mated by the orbital period itself, which is given for 
a pair of galaxies by 



T ■ 



(38) 



where a is the semi-major axis of the pair and m is 
the mass of one galaxy. 

We can apply the usual nav argument for the rate, 
to determine a disruption criterion of the pair. For the 
pair to be disrupted, we have 



navT > 1, 



(39) 



as an approximate formula. Here we can estimate the 
density n to be of order 



= N/R^. 



(40) 



A typical relative velocity value is given by the virial 
theorem as 



2DNm 
R ' 



(41) 



We neglect gravitational focusing for the perturbers, 
which will fly by with a velocity much higher than the 
relative velocity of the two merger candidates around 
the apocenter of their flrst orbit. This means that we 
can use the geometrical cross section 

a = nr^ (42) 

We can make a quantitative estimate for this cross 
section as follows. The average amount of energy 
transferred by a third galaxy which passes at a dis- 
tance r from one of the galaxies in the pair is given 
in the impulse approximation by 

This energy change must be larger than the binding 
energy of the pair jla. Thus, we obtain 



am/v^ 



(44) 



as the maximum encounter distance to disrupt a pair. 

We can find a second expression for this distance, 
using equation |3^. Using both expressions to elimi- 
nate r, substituting v from equation (^), we find 

a>{2/^^f"' RN-^">D^'\ (45) 

as the condition for disruption. Taking into account 
that the first bound orbit, after the initial encounter, 
is likely to be highly eccentric, we can approximate 
the apocenter distance after the first encounter to be 
rapo — 2a. The criterion for tidal disruption from a 
single nearby galaxy then becomes Vapo > Vt^iocai with 

n^iocai - RN-^/'>D^'\ (46) 

Comparing eqs. ^ and we see that the latter 
critical distance exceeds the former by a factor 



(47) 



since > 1 by definition. Typical values for this 
ratio will be ~ 10 for rich clusters. 

In conclusion, we see that the tidal forces exerted 
by neighboring galaxies are less important than the 
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tidal force exerted by the parent cluster as a whole, in 
its effectiveness to disrupt a galaxy pair while it has 
become bound after its first encounter. 

How large is the critical distance rt = T[aYi{rt.iocaU'rt ^global) 
Tt, global for tidal disruption in real clusters? For a 
cluster of N galaxies, and using the notion introduced 
in eq. (p8[), we find that the half- mass radius of the 
cluster is larger than that of the individual galaxies 
by a factor of ND jx^ . Thus, using eq. ( |37| ) we obtain 



(48) 



for the maximum apocenter distance after the en- 
counter at which tidal disruption can be avoided. 

To sum up, tidal effects will disrupt those mergers 
that have a binding energy less than a critical value 
Et, given by 



Et 



rt {NDf/^ 4 



(49) 



where we have used equation g, with E the internal 
energy of one galaxy (E = —1/4 in our units, see 
equation (^)). 

For typical clusters, Et lies in the range 0.01 ^ Et ^ 0.1. 
For example, for a cluster with N ~ 100, D ~ 3, and 
a; = 2, we have Et ^ 0.1, while for a rich cluster with 
TV = 1000, D = 10, and x = 3, we have Et - 0.02. 

5.3. Sensitivity of cross sections to tidal ef- 
fects 

In order to calculate the cross section as a func- 
tion of critical binding energy for tidal disruption Et , 
we could perform a systematic search as we did for 
the case of Et — 0. Since this would require a re- 
run of many of our simulations, separately for each 
value of Et, it would be preferable to check whether 
an appropriate approximation exists that may allow 
us to derive a merging criterion for non-zero Et to a 
reasonable accuracy. 

The merging criterion is expressed as 



where Eq and Ei are the orbital binding energy of two 
galaxies before and after the first encounter, and AE 
is defined as Ei — Eq. AE is a function of collision 
parameters p and u, or rp and described in 

section 3. 

As an approximation, we try the assumption that 
AE is determined essentially by rp alone, at least 
when the relative velocity is close to the critical veloc- 
ity for merging. Thus, we have in first approximation: 



AE = --pv'^rit^ 



and therefore 



*)■ 



(51) 



(52) 



It is reasonable to expact this assumption to be quite 
good for most cases since Vp is determined mainly 
by the potential energy at minimum separation, and 
depends only weakly on v, the velocity at infinity. 

To test whether this approximation is satisfactory, 
we compare the result of our numerical experiments 
and the theoretical estimate obtained assuming that 
AE is independent of the velocity at infinity. Figure 
13 shows the orbital energy of the relative motion of 
the two galaxies, in units of the kinetic energy asso- 
ciated with the critical velocity Vcru 



2E„ 



^orb 



2 



(53) 



El ^ Eo + AE < Et, 



(50) 



plotted against the initial orbital velocity at infinity, 
in units of the critical velocity. Note that Vcrit is cal- 
culated for the pericenter distance rp derived from 
the impact parameter p through the two-body Kep- 
lerian gravitational focusing approximation. The ap- 
proximation of velocity-independent energy dissipa- 
tion, given by equation ( |5^ ) , is given here as the solid 
curve. 

The agreement between theory and experiment is 
excellent for a wide range of initial relative velocities. 
Therefore, in the following we use the assumption that 
AE is a function of rp only. 
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Figure 14 shows the effects of non-zero Et in terms 
on the differential merger rate w^cr. It is clear that 
the reduction in the merging probability is significant, 
especially for the Hernquist model. The reduction is 
largest for lower velocities, for which a typical binding 
energy of the orbit after the first encounter is lower, 
and therefore the vulnerability to tidal effects larger. 

Figure 15 gives the nondimensional merger rate 
R — R{oo) as a function of Et. The dependence of 
i?oo on Et is similar for different models. For rich 
clusters with Et ^ 0.02, the merger rate is 3/4 '--^ 1/2 
of the value for Et — 0. For a more modest clus- 
ter, with Et ^ 0.08, the merger rate has dropped to 
1/3 - 1/5 of the value for Et ^ 0. 

Theoretically, R should reach zero for some finite 
value of Et. Thus, a fitting formula of the form 



R — R( 



E, 



t.o 



Et 



E, 



t.o 



(54) 



would be an appropriate formula. We found Et^ = 
0.25 and 7 = 3 to give a reasonable fit for both a King 
model {Wc = 7) and a Plummer model. 

6. Applications: Merging in Clusters of Galax- 
ies 

As a straightforward application of the merger 
rates we have determined in the preceding sections, 
let us take a cluster of galaxies with a 3-D spatial 
density n, and a 1-D velocity dispersion Ce for the 
motion of the galaxies in the cluster. These galaxies 
are all considered to be given by identical Hernquist 
models, with half-mass radius rh and internal 1-D ve- 
locity dispersion ai. We can use the asymptotic ex- 
pression Rao for the merger rate, given by Eqs. ( |3^ ) 
and (|3^), with the relation between the virial radius 
r„ and the half mass radius r/j given by Eq. (||). 

For a cluster of galaxies, with a one-dimensional ve- 
locity dispersion of ae km/s, we can express our main 
result in physical units as follows. Under the assump- 
tion that all galaxies have identical halos, with inter- 
nal one-dimensional velocity dispersion of ai km/s, we 



find a relative merger rate (per unit volume of IMpc^ 
per Gyr) of 



Rr 



0.0084 



rh 



O.lMpcy V 100 km/s 
(55) 

where n is the density of galaxies in the galaxy cluster. 
If we approximate the cluster by a collection of N 
galaxies evenly distributed within a radius R, we find 
a total merger rate of 



Rtr.t = 0.0020 TV" 



1 Mpc 
R 



J \0.1 Mpc 



100 km/s 
(56) 

corresponding to a fiducial merger time scale of tm = 
5.0 X lO^^yr. In order for merging to occur at a rate 
that is significant over a Hubble time, we would re- 
quire a rate that is two orders of magnitudes larger. 
For example, we could take a cluster with 100 galax- 
ies, within a radius of 1 Mpc. If we keep the same 
fiducial values for the other parameters, given in Eq. 
(^6|), mergers would occur at a rate of, on average, 
once every 50 Myr. 

This merger rate is an overestimate, however, since 
we have not yet taken into account the disruptive ef- 
fects of tidal interactions, discussed in the previous 
section. The corresponding correction factor, stem- 
ming from a finite Et, is fairly large. For example, for 
N = 100, D — 10 and crijug — 0.33, the correction 
factor i?/i?o is around 0.25. Thus, we conclude that 
in the above case a merger occurs only once every 200 
Myr. 

Here we assumed that the cluster is in dynamical 
equilibrium, which is certainly not true for clusters 
with, for example, substructures. To apply our result 
to such clusters, we need rather detailed knowledge 
about how the substructures evolve. 

7. Summary 

We have presented cross sections (§3) and reac- 
tion rates (§4) for merging events to happen during 
encounters of equal-mass spherical galaxies, first in 
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isolation, and then in the presence of a background 
cluster (§5). Our results are easily applied to esti- 
mate the merger rate in clusters of galaxies (§6). 

At first glance, the restriction to spherical galax- 
ies might seem somewhat restrictive. However, our 
numerical calculations as well as our very successful 
analytic approximations show that the merger rate 
is largely determined by the outer few-to-ten percent 
of the mass distribution in the galaxies. For realistic 
galaxies, even with only a relatively small amount of 
dark matter, it is therefore the halo mass distribution 
that determines the merger rate. Even if halos were 
significantly flattened, it would be rather surprising 
if our results would be affected by more than a fac- 
tor of two. For practical purposes, therefore, we are 
confident that our results can be applied to galaxy 
clusters, directly as they are given in §6. 
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Fig. 1. a) Nuclear reaction rates in the core of a 
star, as the product of cross section a and the rate fac- 
tor vf{v), with V the relative velocities of the nuclei, 
and f{v) their Maxwellian distribution, b) Merger 
rates in a rich cluster of galaxies, similar to the stel- 
lar evolution case. 

Fig. 2. — For each value of the pericenter distance 
Vp, circles indicate the choice of asymptotic veloci- 
ties V for which a galaxy encounter has been car- 
ried out. Through linear interpolation of the energy 
dissipation, the critical velocity v{rp) for producing 
parabolic outgoing orbits is found for each Vp value. 
The full line connecting these points indicates the bor- 
der of the region within which merger takes place. 

Fig. 3. — Total energy E{t) of the two galaxies after 
the first encounter, as a function of time t, as specified 
in Eq. (7). The initial lowering of the total energy 
is caused by the fact that the tidal interaction terms 
have been neglected in the calculation of the poten- 
tial. 

Fig. 4. — Maximum impact parameter Pcrit{v) for 
which merging takes place, as a function of velocity 
at infinity v, for Plummer model initial conditions. 
Fig. 4b is an enlargement of fig. 4a, for small impact 
parameter values. 

Fig. 5. — Maximum velocity WcHf (p) for which merger 
takes place. The impact parameter p and Vcrit are 
both asymptotic values, in the limit of infinite initial 
galaxy separation. Solid, short dashed, long dashed, 
short dot dashed and long dot dashed curves are re- 
sults for Plummer (rcui = 23), King (Wc = 1), King 
(Wc = 7), Plummer ( 

i~cut — 4), and Hernquist mod- 
els, respectively. 

Fig. 6. — As figure 5 but with the merger velocity 



'^crit 

{vp) now plotted against the pericenter distance 

Tp. The tidally impulsive limit, the line Vcrit ^ Tp 
is given here for comparison. 

Fig. 7. — Same as figure 6 but including the theoret- 
ical estimate for large Vp. Fig. 7a shows the exper- 
imental results for a Plummer model, together with 
the theoretical estimate Vcrit oc . Fig. 7b shows 
similar results for a Hernquist model, together with 
the theoretical estimate VcHt oc 

Fig. 8. — Merging cross section cr('(;) where v is the 
initial encounter velocity (at infinity). The various 
lines have the same meaning as those in Fig. 5. 

Fig. 9. — The differential merging rate v^a{v) plotted 
as a function of the velocity v for five different galaxy 
models. The curves have the same meanings as in 

figure 5. 

Fig. 10. — Merger rate R{x) as a function of the 
ratio of velocity dispersions: CTg indicates the velocity 
dispersion of the galaxies inside the galaxy cluster, 
and CTj the internal velocity dispersion for the stars 
inside each galaxy model. 

Fig. 11. — As Fig. 10, with the Plummer model and 
the Hernquist model, together with their respective 
fitting formulas Rp{x) and Rh{x), given in Eqs. (34- 
35). 

Fig. 12. — The cumulative mass is plotted here for 
the Plummer model and the Hernquist model. The 

more extended mass distribution of the latter is re- 
sponsible for the higher merger rate of the Hernquist 
model in Figs. 10 and 11. 

Fig. 13. The normalized orbital energy £0^6 (equa- 
tion (52)) after the first encounter plotted against 
V I Vcrit , the initial relative velocity at infinity in units 
of the critical velocity for merging. Open circles are 
the results for the Plummer model and filled circles 
are those for the Hernquist model. The curve shows 
the approximation that AE is independent of v (equa- 
tion (51)). 

Fig. 14. — Same as figure 9 but for non-zero val- 
ues of the critical binding energy for tidal disruption 
Et . Solid, short-dashed, long-dashed, and dot-dashed 
curves are results for Et = 0, 0.02, 0.04, 0.08, respec- 
tively, (a) Plummer, (b) King model with Wc = 1, 
(c) Hernquist model. 
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Fig. 15. Merger rate R{x) as a function of the 
critical binding energy for tidal disruption Ef. Solid, 
short-dot-dashed, long-dot-dashed, and short-dashed 
curves are results for Plummcr (rcut = 23), King 
{Wc = 1), King {Wc = 7), and Hernquist models, 
respectively. 
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